Search for: All records

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide essential insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials 

We investigate the quantum dynamics of a spin coupling to a bath of independent spins via the dissipaton equation of motion (DEOM) approach. The bath, characterized by a continuous spectral density function, is composed of spins that are independent level systems described by the su(2) Lie algebra, representing an environment with a large magnitude of anharmonicity. Based on the previous work by Suarez and Silbey [J. Chem. Phys. 95, 9115 (1991)] and by Makri [J. Chem. Phys. 111, 6164 (1999)] that the spin bath can be mapped to a Gaussian environment under its linear response limit, we use the time-domain Prony fitting decomposition scheme to the bare–bath time correlation function (TCF) given by the bosonic fluctuation–dissipation theorem to generate the exponential decay basis (or pseudo modes) for DEOM construction. The accuracy and efficiency of this strategy have been explored by a variety of numerical results. We envision that this work provides new insights into extending the hierarchical equations of motion and DEOM approach to certain types of anharmonic environments with arbitrary TCF or spectral density.